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Abstract. We compare the slow dynamics of irreversible gels, colloidal gels, 
glasses and spin glasses by analyzing the behavior of the so called non-linear 
dynamical susceptibility, a quantity usually introduced to quantitatively characterize 
the dynamical heterogeneities. In glasses this quantity typically grows with the time, 
reaches a maximum and then decreases at large time, due to the transient nature 
of dynamical heterogeneities and to the absence of a diverging static correlation 
length. We have recently shown that in irreversible gels the dynamical susceptibility 
is instead an increasing function of the time, as in the case of spin glasses, and 
tends asymptotically to the mean cluster size. On the basis of molecular dynamics 
simulations, we here show that in colloidal gelation where clusters are not permanent, 
at very low temperature and volume fractions, i.e. when the lifetime of the bonds is 
much larger than the structural relaxation time, the non-linear susceptibility has a 
behavior similar to the one of the irreversible gel, followed, at higher volume fractions, 
by a crossover towards the behavior of glass forming liquids. 
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1. Introduction 

Many complex systems, such as glasses, spin glasses, irreversible gels, colloidal gels 
and others, exhibit a complex dynamics, all characterized by a slowing down usually 
leading to a structural arrest. Nevertheless there are significative differences, among 
these systems, often not well clarified. One of the key concept to describe the slow 
dynamics in glassy systems is the concept of dynamical heterogeneities. Here we want 
to classify and compare the above systems by looking at the behavior of the so called 
non-linear dynamical susceptibility, a quantity usually introduced to quantitatively 
characterize the dynamical heterogeneities. In glass forming liquids different definitions 
have been proposed for the non linear dynamical susceptibility [H[2]. Mostly considered 
are the fluctuations of the self and the total intermediate scattering functions. These 
quantities typically have similar behavior, grow with the time, reach a maximum and 
then decrease at large time. This behavior is a consequence of the transient nature of 
dynamical heterogeneities and the absence of a diverging static correlation length. On 
the other hand for systems with quenched interactions such as spin glasses, characterized 
by a diverging static correlation length at spin glass critical temperature, the dynamical 
susceptibility defined as the fluctuations of the time dependent spin-spin autocorrelation 
function is a monotonic function increasing with the time. As the time goes to infinity 
this dynamical susceptibility for a fixed temperature T tends to a plateau whose value 
coincides with the static non-linear susceptibility. Therefore as T approaches the spin 
glass temperature the value of the plateau diverges as the static non-linear susceptibility. 
In irreversible gels the definition of the corresponding dynamic susceptibility is not 
straightforward and one should carefully distinguish different dynamical quantities. We 
have recently shown [3] that in a microscopic model for irreversible gels the dynamical 
susceptibility defined as the fluctuations of the self-intermediate scattering function is a 
monotonic function as in the case of spin glasses and, for each fixed value of the volume 
fraction, its long time limit tends to a plateau whose value coincides with the mean 
cluster size. The value of this plateau therefore diverges at the percolation threshold as 
the mean cluster size. Such finding corresponds to the fact that in irreversible gelation 
the heterogeneities are due to the static nature of the clusters. 

On this basis, we speculate that in colloidal gelation, where clusters are not 
permanent due to the finite bond lifetime, this non-linear susceptibility should show 
a behavior similar to the dynamical susceptibility of the irreversible gel at very low 
temperature and very low volume fraction, where the lifetime of the bonds is much larger 
than the structural relaxation time. At higher volume fractions and temperature, it 
should crossover towards a behavior of glass forming liquids. Here we give some evidence 
based on some molecular dynamics simulations of a model for colloidal gelation [UEJE]. 
Moreover, using this scenario, we interpret previous results found in experimental 
investigations of colloidal suspension [7] and in some molecular dynamics simulations [I]. 
Finally, we show that also in spin glass type of models when the lifetime of the interaction 
is made finite the behavior found is similar to that found in colloidal systems [8]. 
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In the following, first we recall the behaviour of linear dynamical susceptibility of 
glass forming liquids and spin glasses (Sect 12. lj ). Later (Sect J2.2l) we consider the case of 
irreversible gels. In Sect J3. 11 we discuss the case of colloidal gelation and compare with 
a spin glass type of model with annealed interactions (Sec H3.2l) . Finally we analyze the 
emerging scenario and the further developments that this study suggests (SectHJ). 



2. Systems with quenched interactions 

2.1. Spin glasses 

We briefly recall the 3d Ising spin glass where the Hamiltonian of the model is 
(ij) e ijSiSji wit h Si — ±1 Ising spins, and — ±1 quenched and disordered 
interactions. The 3c? model undergoes a transition at a temperature, T$g, with a 
divergence of the static non-linear susceptibility, Xni = jj E^K^^i) 2 ]) where the 
average (• • ■) is over the Boltzmann measure, and the average [• ■ ■] is over the disorder 
configurations. The dynamical non-linear susceptibility was firstly introduced in p-spin 
models, considered as the prototype models of glass formers in mean field pQ, 

X (t) = N[(q(t) 2 ) - (q(t))% (1) 

where q(t) = ^ E; Si(f)Si(t' + 1) and the average (■ • •) is done on the reference time t'. 
In the 3d Ising spin glass, differently from the behavior observed in p-spin models, x{t) 
grows monotonically until a plateau value is reached. The plateau value coincides with 
the static non-linear susceptibility and diverges at the transition [8]. 



2.2. Irreversible gels 

In this section we present a molecular dynamics (MD) study of a microscopic model 
recently introduced [3j for irreversible gels. We consider a 3c? system of N = 1000 
particles interacting via a Lennard- Jones potential, truncated in order to have only the 
repulsive part: 

tjLJ = f 4e[(a/r# 2 - (a/r^ + J], n 3 < 2^a 
ij ~ \ 0, Tij > 2 1/6 cr 

where is the distance between the particles i and j. After a first equilibration, we 
introduce quenched bonds between particles whose relative distance is smaller than Rq 
by adding an attractive potential: 



tjFENE 
ij 



-0.5k R 2 ln[l - (r^/Ro) 2 ], r - < R 

OO, Tij > Rq 



representing a finitely extendable non-linear elastic (FENE) [9]. The system is then 
further thermalized. We have chosen k$ = 30e/a 2 and -Ro — 1-5(7 as i n Ref. [9] 
and performed MD simulations in a box of linear size L (in units of cr) with periodic 
boundary conditions. The equations of motion were solved in the canonical ensemble 
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(with a Nose-Hoover thermostat) using the velocity- Verlet algorithm [TU] with a time 
step At = O.OOlSr, where 5r = o^m/e) 1 / 2 , with m the mass of particle. In our reduced 
units the unit length is a, the unit energy e and the Boltzmann constant ks is set 
equal to 1. The temperature is fixed at T = 2 and the volume fraction = 7ra 3 N/6L 3 
is varied from = 0.02 to — 0.2. By varying the volume fraction we find that 
the system undergoes a percolation transition in the universality class of the random 
percolation. We have investigated the dynamical evolution of the system by studying 
the self intermediate scattering functions, & s (k, t) = -k ^2 i=1 e lk '^ n ^'~ n ^\ for increasing 
volume fractions. As in experiments on real gels, we observe stretched exponential 
decays at volume fraction lower than the percolation threshold C , and approaching C 
the onset of power law decays. 

We have moreover measured the dynamical susceptibility associated to the 
fluctuations of the self intermediate scattering functions [2], i.e. X4,(k,t) = 
N [(\§ s (k,t)\ 2 ) — (§ s (k,t)) 2 ], where (...) is the thermal average for a fixed bond 
configuration and [. . .] is the average over the bond configurations. We have shown 
analytically in Ref. [3] that this quantity, in the thermodynamic limit, for t — > oo and 
k — > 0, tends to the mean cluster size. 

In the Main frame of Fig{TJ x^mim t) (with k min = 2n/L) is plotted for increasing 
values of the volume fractions < C . Differently from the non monotonic behavior 
typically observed in glassy systems, we find that it increases with time until it reaches 
a plateau, whose value increases as a function of 0. 

In the Inset of FigHJ Xas(k min , 0) = lim^oo x^min, t) is plotted as a function of 
(0 C — 0) together with the mean cluster size. We find that, as the percolation threshold 
is approached from below, Xas(k min , 0) diverges as a power law at C . The exponent, 
within the numerical accuracy, is in agreement with the value of the exponent 7 of the 
mean cluster size. This finding confirms that one key difference between irreversible 
gelation due to chemical bonds and supercooled liquids close to the glass transition 
is that in irreversible gelation the heterogeneities have a static nature (clusters). The 
clusters, on the other hand, affect the dynamics and as a consequence the dynamic 
transition coincides with the static transition, characterized by the divergence of a static 
correlation length (linear size of the clusters). 

The behavior observed in the case of permanent gels is very similar to that of spin 
glasses in finite dimensions. Although these systems have a very different structures, 
they show a very similar dynamical behavior due to the static nature of heterogeneities. 
In both cases interactions are quenched. What can we expect in the cases where the 
interactions are not quenched and have a finite lifetime? In the following section we will 
try to answer to this question. 
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Figure 1. (Color online)Main frame: Xi(kmin,t) as a function of t for <fi = 0.02, 
0.05, 0.06, 0.07, 0.08, 0.085, 0.09, 0.095, 0.10 (from bottom to top). Inset: Asymptotic 
values of the susceptibility (full triangles), Xas(k m in,4>) an d mean cluster size (open 
squares) as a function of (<p c — 4>). The data are fitted by the power law (<f> c — </>)~ 7 
with 7 = 1.8 ±0.1. 



3. Systems with finite bond lifetime 

3.1. Colloidal gelation 

In this section we present the results obtained in a MD study of a DLVO-type potential 
[TT] for charged colloidal systems, and discuss how in colloidal gelation the finite 
bond lifetime affects the dynamics and, in particular, the behavior of the dynamical 
susceptibility. 

We consider a system of N = 10000 particles, interacting via a DLVO-type 
potential, which contains a Van der Waals type interaction plus an effective repulsion 
due to the presence of charges: 



where a\ = 2.3, a 2 = 6, a 3 = 3.5, and A = 2.5. With these parameters the repulsion 
term dominates the Van der Waals attraction at long range, providing a short range 
attraction and a long range repulsive barrier. The potential is truncated and shifted at 
a distance of 3.5er. To mimic the colloidal dynamics, we performed MD simulations at 
constant temperature. Equations of motion were solved in the canonical ensemble (with 
a Nose-Hoover thermostat) using a velocity Verlet algorithm with a time step of O.OOlio 



(where to — y ^t - and m is the mass of the particles). 

We find jU O [6] that at low volume fraction compact stable clusters form with 
typical size s ~ 10. By increasing the volume fraction a residual attractive interaction 




(2) 
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Figure 2. (Color online) The dynamic susceptibility, Xi(kmin >t) , f° r ksT = 0.15e 
and cj) = 0.10, 0.11, 0.12, 0.13 (from left to right). 

between the clusters produces elongated structures, which finally order in a columnar 
phase [5]. A small degree of polidispersity is introduced [6] in order to avoid the 
transition to the ordered phase. In this case, by increasing the volume fraction the 
elongated structures form instead a long-living random percolating network, i.e. the gel 
phase. The bond lifetime has a non monotonic behaviour: At /c#T = 0.15e it decreases 
of about one order of magnitude from <p = 0.10 to 0.13, has a minimum at 0.13, and 
finally increases for > 0.13. 

The dynamical susceptibility, Xi(^min,t), is measured for = 0.10, 0.11, 0.12, 0.13 
(see Fig. [2]). The data give evidence of a clear crossover from the low volume fraction 
regime to the intermediate regime. In the first regime, where the bond lifetime is much 
larger than the structural relaxation time, we find a behavior resembling that observed 
in the permanent gels: Although Xii^mimt) is a non monotonic function, it increases 
until a value comparable to the mean cluster size is reached; A plateau decreasing slowly 
as a function of time is clearly present in the intermediate time region; Finally, at very 
long times, X^mini t) decreases to its equilibrium value. Increasing 0, the bond lifetime 
and the structural relaxation time become comparable: In this case X^minit) i which 
is again a non monotonic function, displays a well pronounced maximum as usually 
observed in glassy systems. This data suggest that in the first regime the clusters behave 
dynamically as made of permanent bonds as in chemical gelation, and the dynamics is 
dominated by the presence of such clusters. Increasing 0, the structural relaxation 
begins to be affected also by the crowding of the particles, and a clear crossover to a 
new glassy regime is found. 
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Figure 3. Dynamical non-linear susceptibility in the annealed FLG for densities 
p = 0.52, 0.53, 0.54, 0.55, 0.56, 0.57, 0.58, 0.59, 0.60, 0.61 (from left to right). 



3.2. A spin glass type of model with annealed interactions 

It is now extremely interesting to analyse the case of annealed interactions in a spin glass 
type of model. To this aim we consider the results obtained in Monte Carlo simulations of 
the so called frustrated lattice gas (FLG). This model, recently introduced in connection 
with the glass transition [T2J [13] , has mean field properties closely related to those of 
p-spin models. Being constituted by diffusing particles, it is suited to study quantities 
like the diffusion coefficient, or the density autocorrelation functions, that are usually 
important in the study of liquids. The Hamiltonian of the model is: 

- (3H = J ^((, r S',.S' ; - l)niTij + fMy^Uj, (3) 
(w) * 
where Si = ±1 are Ising spins, Hi — 0,1 are occupation variables, and = ±1. In 
the case where are quenched variables randomly distributed the 3c? model undergoes, 
a transition of the type of 3c? spin glasses |14] . Here we considered the case where 
the interactions evolve in time, i.e. they are annealed variables [8]. In this case 
dynamical properties strongly resembling those of glass formers and well fitted by the 
mode coupling theory for supercooled liquids are found. 
The dynamical non-linear susceptibility is defined by 

X (t) = N[(q(tf) - (q{t))% (4) 

where q{t) = i J2iSi{t') n i(t')Si(t' + t) n i(t' + t) and the average (■ • •) is done on the 
reference time t'. In Fig. [3], x{t) is plotted for increasing values of the density. The 
same behavior of the p-spin model in mean field [1] and of MD simulations of the 
Lennard- Jones binary mixture [15] is observed: x{t) shows a maximum, x(^*)i ^ na ^ 
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seems to diverge together with the time of the maximum t*, when the density grows. 
For the highest density, the maximum of x(t) decreases, possibly due to the transition 
to an unfrustrated state. Comparing the behavior found in this case with that shown 
in the previous section, we note that here the first regime with a clear plateau in the 
susceptibility is not present. This is probably due to the fact that the interaction 
relaxation time is in this case always comparable to the structural relaxation time. 

4. Conclusions 

By means of the dynamical susceptibility, we have here analyzed the presence of 
dynamical heterogeneities in systems with quenched and annealed interactions. In the 
case of quenched interactions, as it happens in spin glasses, the dynamical susceptibility 
grows monotonically in time until a plateau value is reached. The plateau value coincides 
with the static non-linear susceptibility and diverges at the transition. This behavior 
is in fact also observed in a microscopic model for irreversible gels, where the plateau 
value of the dynamic susceptibility diverges at the percolation transition as the mean 
cluster size. These results confirm that in irreversible gelation the heterogeneities have 
a static nature (clusters). These clusters, on the other hand, affect the dynamics and as 
a consequence the dynamic transition coincides with the static transition, characterized 
by the divergence of a static correlation length (linear size of the clusters). 

With annealed interactions instead, in the case of spin glass type of models, one 
recovers the non monotonic behavior of the dynamical susceptibility, which is typically 
observed in glasses. This is due to the transient nature of dynamical heterogeneities and 
the absence of a diverging static correlation length. We analyze moreover the case of 
colloidal gelation, where the clusters are not permanent due to the finite lifetime of the 
bonds. We find that the dynamical susceptibility is again a non monotonic function, 
and displays at high volume fraction a well pronounced maximum as usually observed 
in glassy systems. Remarkably, at very low temperature and very low volume fraction, 
where the lifetime of the bonds is much larger than the structural relaxation time, the 
dynamical susceptibility shows a behavior similar to the dynamical susceptibility of the 
irreversible gel with a crossover, at higher volume fractions, towards a behavior typical 
of glass forming liquids. These results suggest that in the first regime the dynamics 
is dominated by clusters, made of bonds which can be considered as permanent in 
this time window. Increasing 0, when these two time scales become comparable, the 
structural relaxation begins to be affected also by the crowding of the particles, and a 
clear crossover to a new glassy regime is found. 
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